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1. Introduction 

The theory of noise in quantum transport in mesoscopic systems is a very active field 
of research [H |2]- In addition to the first few moments of the transmitted charge 
the full probability distribution can be adressed, called full counting statistics (FCS). 
The systems usually studied consist of a finite "dot" -region connected to N leads 
which initially are separated from the dot region and have different chemical potentials 
[31 m El E]. After connecting the subsystems the time evolution of the particle transfer 
between the leads is studied. 

In this paper we focus on systems with two quasi-one-dimensional noninteracting 
leads. The "left" lead consists of Mi tranverse channels and the initial state is described 
by a grand canonical ensemble with chemical potential and temperature = PL/ks- 
This lead is connected via a finite dot region to the "right" lead with Mr = M — Ml 
tranverse channels, chemical potential ^ir < and temperature Tr = Pn/ks. 

For noninteracting electrons the calculation of the characteristic function of the 
probability function of the transferred charge can be exactly reduced to the evaluation a 
time dependent determinant over the full one-particle Hilbert space [7j. For the lattice 
systems studied in this paper the dimension of this space is finite before taking the 
thermodynamic limit. Exact results can be obtained for very large but finite systems 
by numerically calculating this determinant. For times t smaller than the time it takes 
the "charge fronts" which move into the subsystems after connecting them to return 
to the connection point after the reflection at the boundaries the results are almost 
independent of the size of the system |8]. 

The long time limit was investigated by Levitov and Lesovik [3] for infinite leads. 
They presented the leading time order result for the logarithm of the characteristic 
function (linear in t) as an energy integral over the logarithm of a determinant of an 
M X M matrix which involves the scattering matrix s(e) for a single particle. For two 
single channel leads, i.e. a strictly one-dimensional system, the 2x2 determinant was 
explicitely evaluated [3j . The result involving the transmission probability T(e) and the 
Fermi functions of the two leads is usually called "Levitov-Lesovik formula" . In order to 
contrast it from the general leading long time approximation involving the determinant 
of a M X M matrix it is called "two single channel leads Levitov-Lesovik formula" in 
the following. For perfect transmission this formula leads at zero temperature to a delta 
function for the probability distribution of the transferred charge incorrectly signalling 
"zero shot noise". It was shown previously that in the case of perfect transmission at 
zero temperature the logarithm of the characteristic function increases logarithmically 
with time leading to a probability distribution of finite width also at zero temperature 

PE]. 

In addition to exact numerical calculations for M = 2 an explicit result for the 
temperature dependent Levitov-Lesovik M x M determinant is presented in terms of 
the eigenvalues of a temperature independent matrix. The result is a sum of terms in 
the form appearing in the two single channel lead Levitov-Lesovik formula and might 
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have guessed from the well known results for the first two cumulants for this general 
case [T]. 

For the strictly one-dimensional case Ml = Mr = 1 exact numerical results were 
obtained for large but finite lattices at zero temperature |8j. The first step to obtain 
the probability distribution of the number of electrons transmitted to the right lead 
was to calculate the time dependence of the one-particle projection operator PR{t) onto 
the right lead. At zero temperature only Pnit) = noPji{t)no enters, where Uq is the 
projection operator onto the initially occupied one-particle states. The eigenvalues 
Pm{t) of -PR,(t) determine the probability distribution wr of the number of particle 
transferred to the right lead. The time dependent entanglement entropy after connecting 
the subsystems can also be simply expressed in terms of these eigenvalues [10] . 

There are approximately Nt = t{fiL — fiji)/{2Ti) eigenvalues Pm(t) ~ T(/ii^-f-27r(m — 
1/2) /t) different from zero and one in the long time limit, where T(e) is the transmission 
probability [8]. The transition region of the finite eigenvalues to the zero eigenvalues is 
not captured by this expression. In order to obtain analytical approximations for the 
exact numerical eigenvalues in this regime the logarithmic correction in the large time 
limit has to be known. 

For finite temperatures the probability distribution is determined by Pnit) in the 
full one-particle Hilbert space. A clever rewriting of the determinantal expression for 
the characteristic function [11] allows a simple generalizion of the zero temperature 
numerical procedure to obtain finite temperature results. In section 2 exact results 
are presented for a strictly one-dimensional tight binding model, i.e. Ml = Mr = 1. 
In the long time limit an accurate analytical approximation for the eigenvalues of the 
temperature dependent operator X{t) introduced by Adamov and Ibanov (AI) [iT] is 
presented which replaces Pnit) at finite temperatures. 

The analytical expression for the eigenvalues of X{t) is extendend to arbitrary 
values of Ml and Mr in section 4 by using a similar rewriting as used by Adamov 
and Ivanov [TT] for the leading order in t result for the logarithm of the characteristic 
function qr corresponding to wr. This derivation also provides a simple derivation of 
the explicit generalized Levitov-Lesovik formula already mentioned. For the special case 
Mr = 2 a comparison is made with analytical results derived earlier [12]. 

In section 5 the generalized long time approximation is elucidated for a model with 
leads which are stripes of equal width and apart from a single site impurity a perfect 
transition region. 

2. Counting statistics for noninteracting electrons 
2.1. General formulation 

In the following we consider a system which consists of a finite "dot" -region described 
by the Hamiltonian Hq"^ connected to the left and right lead with Hamiltonians Ho a 
with a = L,R. The leads are initially separated from the dot region . The number of 
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electrons in the initial state are Ai"g°* and Ng^a- We assume the intial state |^'(0)) to be 
an eigenstate of Hq°^ and the ifo,a 

l*(o)) = i<"')®|i?^'^)®K°'-). (1) 

The time evolution for times greater than zero is described by the Hamiltonian 

H = + Ho,L + Ho,R + Y,ya = Ho + V . (2) 

a 

The term V which couples the leads with the dot region will be specified later. The 
probability distribution that Q electrons are transferred to the right system after time 
t is given by 

WRit, Q) = {mm - - No,R)]m)) (3) 
1 

~ 2^ 

Here Mr ist the particle number operator of the right lead and gR{t, A) is the 
characteristic function. With the particle number operators Mait) in the Heisenberg 
picture it is given by 

g^{t, A) = (<l'(0)|e'^^«(*)e-'^«|$(0)) . (4) 

The assumption that the initial state is an eigenstate of the particle number operators 
was used. For initially grand canonical subensembles with different temperatures and 
chemical potentials 

(a) /K\ 

and Po°* of the same type, which corresponds to a total statistical operator po of the 
generalized canonical form po = g~^°/Zo the averaging yields for the characteristic 
function 

^?^(t,A) = (e^^««e-'^-> , (6) 

where (...) denotes the averaging with the statistical operator po. This result is also 
valid for interacting fermions. 



2.2. Noninteracting fermions 

For noninteracting fermions the characteristic function can be expressed as a 
determinant in the one particle Hilbert space using Klich's trace formula[71 [8] 

gRit, A) = det [1 + (e^A^^We-^^- - l) no] (7) 
= det [1 + a(t)nQ] = det [1 + fiQaij:)] , 

where no = (e'*" + 1)^^ is the Fermi operator. It is determined by the Fermi functions 
describing the initial state. The equality in the second line holds because the inverse of 
no exists, in contrast to the zero temperature case where no is a projection operator. 
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Using e*'^^^'^*) = 1 + (e*'*' — l)Pfi(t) and the definition d{\) = e^^ — 1 the operator 
1 + n^a can be written in the form proposed by AI [TT| 

= {1 + d{X) [(1 - no)PR + noPi?(t)]} e-*^^« 

= [l + d(A)X(t)]e-^^^« . (8) 

— 1/2 ~ 

As Uq exists and no and Pr commute it is more convenient to work with X{t) = 
Uq ^^'^X{t)nl^'^ i.e. 

X{t) = (1 - no)PR + ny^PR{t)ny^ . (9) 
This yields for the characteristic function for arbitrarily large but finite systems 
g^it, A) = e~^^^«det [l + d{\)X{t) 

N„ 
m=l 

Nh 

= 5^c„(t)e^("-^«)^ , (10) 

ri=0 

where is the dimension of the total one-particle Hilbert space and of the one of 
the right lead. They are both finite for finite lattice systems. The coefficients Cm{t) can 
be obtained recursively from the eigenvalues of the AI one-particle operator X{t) 

defined in Eq. ([9]) as described shortly in the appendix. 

Apart from the replacements Nq^r — > Nr, A^tot — ^ Nh and Pmif) Xmit) this 
finite temperature result has the same form as the T = approach which was used as 
the starting point for the exact numerical calculation of the PCS [8J. The probability 
distribution at finite temperatures is given by 

Nh 

WR{t, Q) = Y1 - - Nr)) ■ (11) 

n=l 

To obtain exact results for the FCS one first has to calculate X(t) using the result 
for PR{t) and then obtain its eigenvalues Xm(t). This is done for the simplest case 
Ml = Mr = 1 in the following section. 



3. Exact results for Ml = Mr = 1 

3.1. The model 

In this section we present exact numerical results for the probability distribution w{t., Q) 
for a one dimensional tight binding model with a one site (noninteracting) dot. The 
unperturbed one-particle Hamiltonians of the subsystems are given by 

Na-l 

ho,a = - t\\ ^ {\am){a{m + 1)| + H.c), 

m=l 

hr = Vo\0){0\ (12) 
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The number of sites in the leads are given by A^^^. In \am) the label a takes the value 1 
for a = R and —1 for a = L. The hopping matrix elements in the leads t\\ are taken as 
unity in the numerical calculations which leads to total bandwidth of 4. The eigenstates 
of the unconnected leads are standing waves. The coupling between the subsystems is 
described by the hopping term 

v = -tL\-l){0\-tR\0){l\+H.c. (13) 
3.2. Numerical results 

The first step to calculate the eigenvalues Xm{t) is to obtain PR{t) using the time 
dependence of the one-particle states \ea^) 

Nr 

{e^^^\PR{t)\4^) = • (14) 

m=l 

The time dependence of the states is calculated using the spectral decomposition of the 
full one-particle Hamiltonian [8]. From the resulting A'^^^ x A'^^^ matrix one obtains X{t) 
as prescribed in Eq. (Q. In the following we show results for identical temperatures in 
the initial subsystems and fiL = yUdot > A^i?,- 

We begin with a generic example = 0.8, = 0.5 and Vq = 0.4. As the coupling of 
the dot to the leads is asymmetric the transmission probability at the resonance energy 
is less then one. The results shown are for time t = 200. If twice the number of lead 
sites Na is larger than fmax^ = 2t the results for the eigenvalues X„ which differ from 
zero and one become independent of the Na. For the exact numerical results shown we 
used Nl = Nji = 500. As X{0) = Pr holds there are Nr eigenvalues one and A'^ + 1 
eigenvalues zero at the initial time. Therefore we show the eigenvalues in descending 
order as a function of u — Nr. In Fig. 1 the eigenvalues X„ are shown for three different 
temperatures. 

As mentioned in the introduction the "zero temperature" result (3 = 2000 can be 
well described analytically except for the narrow transition region to the zero eigenvalues 
which is related to the logarithmic corrections in the long time limit [HI [9]. For /3 = 10 
the main effect is to smooth out this transition region. A new effect sets in at larger 
temperatures. For {3 = 2 part of the eigenvalues with n < Nr which are one at the 
initial time get visibly reduced. 

In order to understand this behaviour analytically we start from the Levitov-Lesovik 
formula [3, 4J 

In gR{t, A) = ^ + A))cie , (15) 

where B = 2 for the choice t\\ = 1 and 

F{e, A) = T(e) {d{X) h{e) fR{e) + rf*(A)/«(e)^(e)) , (16) 
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Figure 1. Eigenvalues of X{t) which differ from one and zero for tL = 0.8, tu = 0.5 
and Vq = 0.4 at time t — 200 for three different tempcratm'es, i.e. values of p. The full 
circles correspond to (3 = 2000 which hardly differs from the zero temperature result. 
The open squares are the results for (3—10 and the open triangles for P — 2. 



whith fa = ^ — fa- We approximate the integral in Eq. f[T^ by a finite Riemann sum 
over N intervals of size 2B / N and use the trapecoidal rule 



t 2B 



N 



1 + F -5 + (j 



1 25 



X 



which agrees with the integral in the limit N —>■ oo. As Eq. 
approximation for the large time limit we choose = N{t) with 

t 



(17) 

ffTSl) is itself an 



Nit) 



27r 



-2B 



Then the prefactor in the sum equals one and gR{t) takes a form which can easily be 
compared with Eq. lfTOj) 



N(t) 

^ij('^,A)^n[i+^(^^'^)] ' (19) 

where ej = -B + 27r(j - l/2)/t. 

In the low temperature regime ksTa ^ I^l — IJ^r the approximation for the 
eigenvalues Xj can easily be read off as the factor multiplying d* in F(e, A) in 

Eq. (fT6ll is exponentially small. With e~*'^(l + d{X)) = 1 the comparison with Eq. liTOl) 
shows that the eigenvalues different from one are given by 

X, ^ T(6,)^(e,)/«(e,) ^ r(6,)(/i(6,) - fnie,)) (20) 

At zero temperature this agrees with the result mentionend in the introduction [8j. 

At arbitrary temperatures one has to factor 1 + F{e, A) with F defined in Eq. (fT6ll 
in the form 

1 + F(e, A) = e-'\l + a+(e)rf(A))(l + a_(e)rf(A)) (21) 
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Figure 2. Long time approximation Eq. (|22|) for the eigenvalues X{t) plotted as 
a function of the ej for the parameter values of Fig.l . Filled dots: f3 = 2000, open 
squares: /3 = 10, open triangles: (3 = 2. 



The comparison with Eq. ( 1161) yields 



1 + T{fL - fR) , . 

a± = ± w (22) 



with 



i.^^nfLjn^ -T/^/i? (23) 




TUl - Ir) 



TUl + Ir] 



2 

+ T{l-T)fLfR. 

The second form for w is useful for the discussion of the low temperature results and 
the third form shows that there is an energy gap in the spectrum for non-perfect 
transmission. It also shows that the factorization is simplest for perfect transmission 

1 + F(e, A) = e~'\l + fL{t)d{X)){l + fR{t)d{\)) . (24) 



In Fig. 2 we show the analytical approximation Eq. (122|) for the eigenvalues for 
the parameter values used in Fig. 1 as a function of the ej. The result for f] = 2000 
(filled dots) is almost identical to the zero temperature results. The eigenvalues a_(ej) 
are nonzero in the energy range between /i/j and fi^ and show the energy dependence of 
the transmission probability. The eigenvalues a+(ej) are almost identical to one. Also 
for (3 = 10 (open squares) the a+(ei) equal one within the drawing accuracy while the 
a_(ei) are very well approximated by Eq. (|20l) and show the thermal broadening of the 
zero temperature result. For {3 = 2 (open triangles) Eq. (|20i) no longer presents a good 
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Figure 3. Comparison of tiie exact eigenvalues of Fig.l (shown as the full lines) with 
the approximations a_ presented in descending order. The results for /3 = 10 and 2 
are shown in a restricted range. 



approximation for the a_(ej) and the are a+(ej) are clearly smaller than one in a large 
energy range. 

In order to compare the approximate eigenvalues a± with the exact numerical 
eigenvalues the a^. have to be brought into descending order. This is shown in 
Fig. 3. For (3 = 2 the a_ (open triangles) agree very well with the exact results apart 
from the fact that there are pairs of almost equal eigenvalues. As shown later this has a 
rather small effect for the calculation of the probability distribution wr. For (3 = 10 the 
deviations are a bit larger. As shown earlier the transition region to the zero eigenvalues 
for P = 2000 is not captured by the approximation T(ei)[8]. 

In order to obtain the probability distribution Wpi{t, Q) form the eigenvalues Xn 
the product in Eq. ffTOj) has to be evaluated recursively. The details are presented in 
the appendix. The results for the parameters of Fig. 1 are shown in Fig. 4. For (3 = 2 
the exact result (full squares) agrees very well with the one using the approximation of 
Eq. fl22l) (open squares). For /? = 10 the agreement of the exact result (full circles) with 
the approximation (open squares) is still rather good. For [3 = 2000 the filled and open 
symbols no longer overlap. 

The discrepancy between the exact results and the approximation Eq. ( l22i) is most 
prominent for perfect transmission in the zero temperature limit. While the Levitov- 
Lesovik formula predicts a single delta function of weight one ("zero shot noise") the 
exact result clearly has a finite width [8] as shown in Fig. 5 where the ti = 1, = 1 
and Vo = was used leading to perfect transmission. 
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Figure 4. Weights of the probabihty distribution wii{t, Q) for the parameters of Fig. 
1. The exact results are presented by the full symbols, the approximation using Eq. 

by the open ones: (3 = 2: filled squares, open circles, /? = 10: filled circles, open 
squres, /3 = 2000: filled triangles, open triangles. 




Figure 5. Weights of the probability distribution wji{t,Q) for the case of perfect 
transmission. The exact results are presented by the full symbols, the approximation 
using Eq. ((22|) by the open ones: (3 = 2. filled circles, open squares, (3 = 10: filled 
squares, open circles, (3 = 2000: filled triangles, open triangles. 



3.3. Perfect transmission 

As the discrepancy for very low temperatures concerns the width of the approximate 
distribution shown in Fig. 5 it is useful to dicuss the behaviour of the second order 
cumulant 

K2{t) = (25) 

n=l 

N{t) 
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Figure 6. Exact results (filled symbols) for second order cumulant K2 for perfect 
transmission as a function of time for three different temperatures: (3 — 2000: filled 
squares, /3 = 100: filled triangles and /3 = 10: filled circles. The dotted line shows the 
linear part t/(7r/3) for /3 = 10. 



Nit) 



J2 [TMlIl + fRfn) + T{1 - T)ifJn + Ml)] 



The (leading time order) shot noise contribution [13] proportional to T(l — T) vanishes 
for perfect transmission and the remaining term is well known in the limit where the 
sum is replaced again by the integral [H El H] • 



Because of the factorization in Eq. fl24l) of 1 + F for perfect transmission the 
derivative of Ingpi with respect to iX takes the simple form 



where the factor a in ajja again takes the value 1 for a = R and —1 for a = L. In the 
wide band limit Pa{B — \fj.a\) ^ 1 the logarithmic terms can be neglected and wr is a 
Gaussian in the Levitov-Lesovik approximation with a temperature independent mean 
value [1] . The fact that the average charge transfer for /5 = 2 in Fig. 5 is slightly less 
than for f3 = 10 and /3 = 2000 is related to the correction term in Eq. fl26l) which adds a 
constant contribution for A = 0. 

In the wide band limit In gji{t, A) can be calculated analytically also for an energy 
independent transmission probability which differs from one [14] . 

Exact numerical results for K2 are shown in Fig. 6 as a function of time for system 
size Nl = Nji = 500 up to times were the result is independent of this system size. 
For (3 = 10 (filled circles) there is an almost linear increase of K2 rather quickly. The 




(26) 
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dotted line shows the hnear increase which follows from the finite temperature Levitov- 
Lesovik formula (see Eq. fl26l) ). For /3 = 100 (filled triangles) there is crossover from 
a logarithmic increase to a linear time dependence around t ^ 50. For (3 = 2000 the 
times shown are too small to see the corresponding crossover and the shown increase 
is logarithmic in time. This logarithmic behaviour [8[ IHl 15] is not captured by the 
approximate eigenvalues a± in Eq. fl22|) . 



4. Long time behaviour for general values of Ml and Mr 

To study the long time limit one first takes the thermodynamic limit in which the 
energies of the unperturbed eigenstates \e,a,i) with i = 1, ...,Ma of the left and right 
lead form M^ and Mr continua extending from ej^-j^ to e^!^^ which usually partially 
overlap. For a given energy e the "open scattering channels" in the leads are those 
for which e G [Cmin' ^max]- Their number is denoted as Ma{e) < Ma- For a given e we 
order the channels with i = 1, Ma{e) as the open channels and i = Ma{e) + 1, Ma 
the closed channels. The dimension of the scattering matrix s(e) for a single particle 
[16] is given by M(e) = Mi(e) +MR{e) < M, with Ma{e) > 1 in order to have scattering 
at all. 

In the following we use a Dirac notation in the M(e)-dimensional space with the 
orthonormal basis \a, i) where the i run from 1 to Ma(e). Then the projection operators 
Pa(e) on the lead channels read 

Pa{e)= \a,i){a,i\ . (27) 
1=1 

The energy dependence of the -Pa(e) is only via the selection of the open channels. 

For the general geometry discussed in section 2 the long time linear in t contibution 
to In gji takes the form[3l [T2] 

Ingnit, A) = — / lndet[l + c(e, X)]de . (28) 



27r 

" "-min 

Here c(e, A) is the M(e) x M(e) matrix 

c(e, A) = P(e)e(A, e)s{e)e\X, e) - l] /(e) , (29) 

with 

e(A, e) = PL{e) + e'^Pnie) = 1 + d{X)Pn{e) 

/(e) = f^{e)PL{e) + h{e)Pn{e) . (30) 
The values of emin and emax depend on the details of the microscopic model. 

Now one can essentially repeat the steps used to derive the Al-form Eq. (fTOj) to 
obtain 

det[l + c(e, A)] = e-^^*^«(^Met [l + d{X)X{e)] (31) 
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with 

X = (1 - f)Pn + f/'s^pRsf'/' . (32) 

The determinant in Eq. (|3T1) can be calculated by first solving the eigenvalue problem 
for X(e). As the M(e) x M(e)-matrix X(e) in Eq. (l32|) is temperature dependent it 
looks as if one has to solve a different eigenvalue problem for each temperature. In 
the following we show that this is not necessary. In fact it is sufficient to solve a single 
Mii{e) X M/j(e) eigenvalue problem to obtain a new generalized long time approximation 
for arbitrary temperatures. To show this we write 

X{^)= E [fR\R.3){R.3\ + h){h\] (33) 

with 

|6,) = /'/V|i?,j) , (6,| = (/^,J>/l/^ (34) 
Multiplying the eigenvalue problem X\Xa) = Xa\Xa) from the left with (6j| yields 

J2 [ih\R,j)iR,j\X^)fR + ih\b,)ib,\X^)] =XMXa) (35) 
i=i 

and multiplying with {R,i\ gives for the overlaps {R,i\Xa) 

{R,i\X^) = ^ V {R,i\h,){hj\X^) . (36) 

Inserting this into Eq. fl35l) leads after multiplication with Xa — Jr to 



J2 [fRibi\PR\bj) + {X^-fR){biK 



bj\X^] 



-2 



= {Xi-fRX^){h\X^) . (37) 

This equation can be rewritten using 

{bi\bj) = fRSij + Ul - fR){R,i\spLS^R,j) . (38) 

Therefore a single hermitian operator in the Mij(e)-dimensional subspace spanned by 
the \i) = \R,i) determines the original eigenvalue problem 

[ih - fR)X^ - hh] r|Ff = {Xl - XJirj^)), (39) 

where = pRsf^/^lX^) and 

r = PrsPls^Pr = {PrsPl){PrsPl)^ = AA^ . (40) 

The elements of the matrix A are the left to right transmission amplitudes [T7] . 
The matrix elements of r are given by G [1, M/j(e)]) 

AI{e) 

E ^^Ke)4(e). (41) 

l=MRie)+l 
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After solving the eigenvalue problem T|r^) = t^|t^) the determination of the in 
Eq. fl39l) . after multiplying with (r^|, is reduced to solving a quadratic equation. With 
a ^ fi,± the solution reads 



X^,,^l±MI±^±^(l±I^^±^'^ -rjj,. (42) 

This is like Eqs. ( I22l) and (!23l) with T(ej) replaced by T^(ei). 

For Male) = 1 there exists only one eigenvalue ri, which for Mi(e) = 1 is given by 
|'5i2(e)P = T{e). For MR(e) = 1 and arbitrary values of Mi(e) the single eigenvalue is 
given 

n(e) = X]|.n(e)r. (43) 

1=2 

This corresponds to the simplest generalization of the Levitov-Lesovik formula[8]. 
For Mji{e) = 2 the two eigenvalues of r(e) are given by 




ri,2 = — ±W — +detr. (44) 



Zero temperature results for this case were presented earlier [12] using a different deriva- 
tion. For the special case Mi(e) = 1 the determinant of r vanishes and only one eigen- 
value of T is different from zero. Generally the number of eigenvalues of r which differ 
from zero is less or equal than M<(e), where M<(e) is the smaller of the two Ma{e). 
This stems from the fact that AA'^ and A'^A have the same nonvanishing eigenvalues 
with the same multiphcities. For Mi^e) < MR^e) one better calculates the eigenvalues 
of A which is a hermitian Mi(e) x Mi(e) matrix. 

In order to compare our general result with the Levitov-Lesovik formula one can 

use 

e-'^(l + rf(A)X^,+)(l + rf(A)X^,_) 

= e-^^ [1 + (1 + (/l - fR)r,)dW + TjnfLdiXf] 
= 1 + T,idiX)fjR + diXrfnh) ■ (45) 
With det(l + d{X)X^{e)) = UaC^ + d{X)X^{e)) this yields 



Ingn^— In [l + r,{d fjn + d*fRh)] de. (46) 

This completes the derivation of the generalization of the Levitov-Lesovik formula for 
two general quasi-one- dimensional leads [H] . The integration has to be split up into iV/ 
energy intervals from e*^*-* to 6*^*+^-' where i = 1, Ni and e'-"^-' = emin and e*^^-^"*"^) = emax- 
In the intervals the number of open channels determined by the Ma(e) is constant. We 
denote the constant value of M<(e) in the m-th interval by M^^\ The splitting is 
discussed in the following section for simple model system. 
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With the eigenvalues X„(e) of the M{e) x M(e) matrix X{e) and the transition from 
the integral to a finite Riemann sum as in Eq. f|T7|) one can obtain an approximation for 
the eigenvalues of the operator X{t) defined in Eq. (j9]) which generalizes the introduction 
of the a± in section III. The approximation for g^i reads 



9RitA)-^ll^-'''''<^'''-^'^ U n [l + rf(A)X.(6f))] (47) 

with Nm{t) = (e*^'"^"'^^ — e*^'^^)t/27r and the energy variables are given by e^™"* = 
g(m) _|_ 27r(j — 1/2) /t. Therefore the approximation for the eigenvalues Xjn{t) which 
differ from 1 and are given by the Xa{e^™'^). The time dependence enters via the ej™"*. 

5. Almost perfect stripe 

In order to elucidate our general result Eq. P6|) we consider leads which are both stripes 
of width N±. Analytical results for the eigenvalues are presented for the case where 
the dot region is identical to the leads except for a single site impurity 

oo 

h = — [ t\\\m,n){m + l,n\ + t±\m,n){m,n + 1\ 

m=—oo n=l 

+ + K)|0, no) (0,no| = /lo + ^o • (48) 

The tildes are introduced to indicate that the separation of the Hamiltonian in the 
unperturbed part and the perturbation is different from the one used in Eq. ([2]) and Eqs. 
^ and (US]). 

For Vo = this is an ideal infinite stripe with eigenvalues 

efc,/ = -2t\\ cos k - 2t^ cos = eJJ + , (49) 

iV_L + i 

with k G [— TT, vr] and / = 1,...,A^. For a given energy e there are 2A^(e) < 2A_l 
scattering channels open for which |e — e^l < 2|t||| holds. The standing wave lead states 
\e,a,l) are labeled hj a = L,R and the transverse quantum numbers / of the open 
channels. The scattering matrix is obtained via the t-operator \161 

t{z) = V + vg{z)v , (50) 

where g{z) = {z — h)^^ is the exact resolvent and v is the generalization of the operator 
connecting the leads with the dot in Eq. (fT3l) . For the almost perfect stripe it is given 
by 

v= -t\\ ^(1 - l,n){0,n\ + \0,n){l,n\ + H.c) 

n=l 

N_L 

= -hJ2i\k-i)){ho)\ + \ho)){hi)\ + H.c) , (51) 
1=1 
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where the |/(m)) are the standing wave eigenstates in the perpendicular direction formed 
from the states \m,n). Then the t- matrix elements take the simple form 

(e, a, l\t{z)\e, a, I') = tj(e, a, ^|/(a))(/(o)|5'(^)l^(o))(^(a')l^' (52) 

where \l(a)) is the standing wave state at n = — 1 for a = L and n = 1 for a = R. The 
t-matrix elements enter the scattering matrix for ^ = e + iO. 

The exact resolvent matrix elements can easily be calculated for the Hamiltonian in 
Eq. fHHl) as the site impurity provides a separable perturbation. With go{z) = (z — /iq)"^ 
one obtains 

(^(o)l^l^'(o)) = (^(o)l^ol^'(o)) (53) 
(/(o)|^o|0, no)Vo(0, nol^ol^(o)) 
1 -Vo{0,no\go\0,no) 

For the open channels 

(/{o)|^o(e + ^0)|r(o)) = 5u' I ~' = (54) 

iBf, - (e - eif 



holds, with B\\ = 2t\^. The gQ-mathx element in the denominator in Eq. (l53l) involves 
contributions from the open and the closed channels 

-i|(/(o)|0,no)P 



(0,no|^o(e + .0)|0,no)= ^^P^^^ (55) 

Z(open) y -^W ~ ~ 

|{'(o)|0,r!„)p 



l{closed) Y ~ 

^vrpoo + 9oo ■ 



?-B^ 



While the contribution of the open channels is purely imaginary the one of the closed 
channels is real. To complete the calculation of the t-matrix elements in Eq. flS^ the 
overlaps (e, a, l\l{a)) are needed. They are related to the density of states at the boundary 
of a semi-infinite chain 

t||(e,a,/|/(,)) = (^y/5;-(e-Q)2) ' (56) 

This leads to the scattering matrix 

Sai,a'u{t) = 6ii'{6aa' " 1) + {l\no)u{no\l') (57) 

with 

{l\no) = = , u = ■ (58) 

/Bf - (e - q)2 1 ~ ^0^00 + ^ttVoPoo 



With the iV_|_(e) x N^{e) projected scattering matrix srl = — 1 + |no)u(rao| the operator 
T defined in Eq. (1401) is given by 

r = 1 - (m + M* - {no\no)\u\'^)\no){no\ . (59) 
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Because of the separable form of r — 1 the only eigenvalue of r different from 1 is given 
by 

Ti = l- {u + u* - {no\no)\u\'^){no\no) . (60) 

Using {no\no) = Z]i(open)(^o|0(^l'^o) = Poo/2 one finally obtains 

^ (1 - Vo& . . 

' (1 - Vo~g§,r + Wpoor ' ^ ^ 

T"2 = T3 = ••• = TAfxCe) = 1 ■ 

The "perfect transmission" eigenvalues Tj = 1 yield contributions to In g^i of the form 
discussed following Eq. fl26|) . As the energy integrations in Eq. (j^6l) for fixed M<(e) are 
over restricted energy ranges the logarithmic corrections in Eq. fl2B]) are important. 

The simplest case is a N± = 2 ladder system. For the special case = t± = 1 there 
are two bands corresponding to e/" = ±1 of width 4. Therefore one has to distinguish the 
energy intervals [—3, —1], [—1, 1] and [1,3] with one, two and one open channel. When 
both channels are open T2 = 1 and g^Q vanishes. This implies ri — >■ for Vq — >■ 00, while 
Ti stays finite in this limit when only one channel is open. 



6. Summary 

In this paper we have generalized the exact numerical method to calculate the FCS for 
large but finite systems [8] to finite temperatures using the eigenvalues of the operator 
X{t) in the Hilbert space of a single particle introduced by Abanov and Ivanov[TT]. 
In the long time limit the results for the probability distribution for the number of 
transmitted particles agree well with the result using the Levitov-Lesovik approximation 
[21 H] except for the case of (almost) perfect transmission. 

Using a similar identity for the finite temperature leading order in time result for 
the logarithm of the characteristic function a new explicit result for In gji was presented 
in Eq. ( l46l) for two general quasi-one- dimensional leads which involves the eigenvalues of 
a matrix formed from the transmission amplitudes. For a simple model these eigenvalues 
were calculated analytically. 
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Appendix 

Recursive step in the calculation of WR{t, Q) 

The numerical finite temperature results presented in section 3 were obtained by first 
calculating the eigenvalues X^it) and then performing the product in Eq. (fTOj) . This is 
done iteratively as follows. 
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Let -F/v(x) be a polynomial given in the form of a product 

N N 
i=l m=0 

The coefficients c^,f ■* are obtained iteratively by calculating the polynomials Fm{x) with 
coefficients ci^'' starting with M = 1 and using 

FM+iix) = (flM+i + bM+ix)FMix) (.2) 

M 

This leads to the recurrence relations 

JM+l) _ AM) ( ^^ 

4*'"^'^ = «M+ie) + &M+i4i\ , 1 < m < M 

^(M+l) _ , (M) 
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